Measuring Condensate Fraction in Superconductors 
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An analysis of off-diagonal long-range order in superconductors shows that the spin-spin correlation 
function is significantly infiuenced by the order if the order parameter is anisotropic on a microscopic 
scale. Thus, magnetic neutron scattering can provide a direct measurement of the condensate 
fraction of a superconductor. It is also argued that recent measurements in high temperature 
superconductors come very close to achieving this goal. 
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OfT-diagonal long-range order (ODLRO) [y^] provides 
an intriguing characterization of superconductors and su- 
perfiuids for which the order parameter signifies a unique 
coherence property of a macroscopic quantum system 
corresponding to the spontaneous breakdown of global 
gauge symmetry. It is therefore unfortunate that direct 
measurements of this order parameter are few and far 
between 

For BCS superconductors, all definitions of supercon- 
ductivity are directly linked to the energy gap in the sin- 
gle particle excitation spectrum. Therefore, the observa- 
tion of the gap, which is abundant, is virtually the same 
as the observation of ODLRO. In contrast, the high tem- 
perature superconductors sometimes exhibit an appar- 
ent gap at temperatures well above the superconducting 
transition temperature, Tc, similar to the gap seen be- 
low Tc In those instances, the presence of a gap is 
not the same as the presence of ODLRO, and it becomes 
necessary to explore ODLRO in more general terms, not 
restricted by the BCS theory. 

The purpose of the present Letter is to offer some in- 
sight into the question of ODLRO in high temperature 
superconductors and to indicate experiments that may 
shed light on this topic. In fact, we argue that recent 
magnetic neutron scattering measurements come tanta- 
lizingly close to accomplishing this goal 1^. We show 
that the instantaneous spin-spin correlation function is 
influenced by the existence of ODLRO. This is partially 
because the spin operators are composite Fermion opera- 
tors, and the correlation function can acquire anomalous 
expectation value, but also because the anisotropic order 
parameter strongly enhances the effect. We also deduce 
that strong intercation effects are necessary for the mag- 
nitude of the effect to be observable. 

a. Summary of experiments: We begin with a brief 
summary of experiments. The magnetic neutron scat- 
tering cross section is proportional to the equilibrium 
dynamic structure factor, 5(k,a;). The corresponding 
instantaneous structure factor, 5(k) = J^^dujS(k,ijj), 
is 



where vq is the volume of the unit cell, and S is the spin 
operator. 

For a class of bilayer cuprate superconductors, 
YBa2Cu306+5 and Bi2Sr2CaCu208+5, corresponding to 
optimal doping, the 5(k, lu) has a resolution limited res- 
onance peak, at an energy « 40 meV and the special 
wavevector k = q = (tt/o, tt/o, tt/c), for T < T^, a and c 
are the lattice constants. More precisely, the 5(k, oj) can 
be classified into parts that are even and odd with re- 
spect to the individual layers of a bilayer. The resonance 
appears in the odd channel (q^ = '^/c), and its energy 
is almost independent of temperature, but its intensity 
vanishes for T > Tc. For underdoped superconductors, 
there is also a peak for T > Tc, but its character is quite 
different, as we comment below. 

The existence of the resonance peak below Tc allows us 
to isolate a distinct part of the instantaneous structure 
factor by integrating with respect to frequency across the 
resonance, which has a simple theoretical interpretation 
in terms of the condensate fraction, as we shall see. In 
this sense, the resonance peak plays an essential role; 
in its absence, it would be difficult to obtain the same 
information, as in La2-iSr5Cu04. 

b. The Off-diagonal long-range order: The ODLRO 
for a singlet superconductor is defined as a property of 
the two particle density matrix, which is 

t, 2 i;3 T,4 i) = (V4(xi)V;I(x2)Vi(x4)^T(x3)) (2) 

where is the creation operator of the Fermion field, 

of spin up and location x. Its existence implies that 

t,2 i;3 T,4 i) - <i>(xi t,X2 i)*<i>(x3 T,X4 i)| 



< (^)7(|xi 



X3I, |X2 - X4I), (3) 



5(k)=t,oEe''"'''(S(0).S(l)), 



(1) 



where the the non- negative function 7(?'i,r'2), which is 
independent of N, tends to zero if either ri or r2 tends 
to 00. 

Following Rcf. 0, one can show that in the thermo- 
dynamic limit, the number of particles N —i- 00 and the 
volume V ^ 00 such that (N/V) n, the function <I> 
is the eigenfunction of the two particle density matrix 
corresponding to the largest eigenvalue Am H 7 given by 
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Am = J dxidx2|$(xi t,X2 (4) 

where Xm/N is the condensate fraction defined to be the 
fraction of electrons participating in pairing, not to be 
mistaken to be the superfluid density, n^, which is the 
stiffness with respect to an imposed twist of the order 
parameter. Hence, in the thermodynamic hmit, the fol- 
lowing spectral decomposition of the density matrix holds 
for all separations: 

t,2 i;3 T,4 i) = <i>(xi t,X2 i)*<&(x3 T,X4 i) 

+ G(xi t,X2 i;x3 t,X4 i). (5) 

The function G vanishes when the separation between 
the groups (xi,X2) and (x3,X4) tends to oo. While the 
first term of Eq. ^ vanishes above Tc, the second term 
continues smoothly to the normal state above Tc. 

For a translationally invariant system corresponding 
to a condensate of zero total momentum, we can write 
$(xi T,X2 i) = <f(xi -X2)(|n)-UT))/V2, where 
$(x) = $(-x). 

For a BCS superconductor (in an abbreviated nota- 
tion), 

p^^) ^ $o(xi2)*$o(x34) + go(xi3)*go(x24), (6) 

where the BCS function <I>o(xi2) = $o(xi — X2) is a func- 
tion of the separation between the mates of a Cooper 
pair. Its integral (1/iV) / dx|$o(x)P is the condensate 
fraction. In BCS theory, and in the absence of disor- 
der, Us — n, while the condensate fraction is ~ iV(0)A, 
where A is the gap, and iV(0) is the density of states at 
the Fermi energy. As T ^ Tc, the condensate fraction 
tends to zero as A'^N{0)/Tc, similar to Us- In general, 
these are distinct concepts, however. 

The BCS function 50 (x) is the parallel spin correla- 
tion function modified by the presence of a gap, which 
smoothly continues to the normal state. In contrast, the 
function <I'o(x) vanishes in the normal state. 

The discussion above can be summarized as follows: 
(1) the spectral decomposition in Eq. (^) is general and 
does not depend on the BCS theory; (2) the order pa- 
rameter function $ includes, in principle, all effects of 
electron-electron interaction; (3) in common with BCS 
theory. Am vanishes above Tc, while the contribution due 
to G continues smoothly above Tc, (4) as in BCS theory, 
the function G must be a global gauge singlet. 

c. The spin-spin correlation function: From here on 
we shall abandon the continuum notation and adopt the 
more natural lattice notation. Assuming that the spin 
rotational symmetry is unbroken, we note that 

(S(li) • S(l2)) = ^<5i,2no - T, 2 i; 2 t, 1 I), (7) 

where Uq is the average occupation of a site. Making use 
of the spectral decomposition of the two-particle density 
matrix, we can write 



(S(ii) • S(i2)> - Ik^.uo - ^[<f>(ll T, I2 i)*<I'(l2 T, ii i) 
+ G(iiT,i2i;i2T,iii)]. (8) 

Then remembering spin rotational symmetry and using 
Eqs. P and |l], we get 

5(k) = ^vono - \v,Y.^'^'^'V^m\]^)? + G(l)], (9) 
1 

where /(I) is now assumed to be normalized 
(woEil/(l)P = !)■ From the perspective of ODLRO, 
the summation variable 1 in Eq. is a degree of freedom 
internal to the Cooper pair, denoting the separation be- 
tween its mates. In the absence of simultaneous diagonal 
long range order, G(l) is a short ranged quantity, which is 
qualitatively unaffected by superconductivity, essentially 
because the short-distance electron-electron interaction 
is similar in both the normal and the superconducting 
phases. The underdoped cuprates, discussed below, are 
more complex. 

We first consider a single CuO-plane, ignoring the fiuc- 
tuations that destroy ODLRO at any finite temperature, 
T; the coupling between the planes is discussed below in 
section (/). The wavevectors k and q are to be inter- 
preted as two-dimensional (2L)) vectors, strictly k|| and 
q||, until the section (/); similarly, Ri is to be interpreted 
as a 2D vector. We hope that this is not a cause for con- 
fusion. 

d. Optimally doped superconductors: Consider first 
the case of optimally doped superconductors for which 
5(k) is peaked at k = q for T < Tc, but not so for 
T > Tc- Here, q is the vector (7r/a,7r/a). In Eq. (j^), the 
exponential factor changes sign rapidly on the scale of a 
lattice spacing when k = q . For any functions 
and G(l) that are smooth on the scale of a lattice spacing, 
the sum will be negligible. 

An order parameter function /(I) such that van- 
ishes when the vector Ri connects sites belonging to the 
same sublattice of a bipartite square lattice [[7| can lead 
to a peak in 5(k), because 

^e"i-^'|/(l)P^-5]|/(l)r, (10) 

1 IGA 

where the sum on the right hand side is over the sublat- 
tice A. Clearly, a d-wave order parameter satisfies this 
criterion, but so does an anisotrpic s-wave order. The 
distinction between them can only be made on energetic 
grounds. The range of /(I) is the "size of the Cooper 
pair" . 

From the argument given above, we expect G(l) to be 
largely unaffected by superconductivity. We can then 
focus on the quantity 

A5(k + q) = f«oAMEe''''''l/(l)r (11) 
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Then, noting that /(I) is normahzed, we get 

ASjq) ^ 3 Am ,^21 
N 2 N ' ^ ' 

The inverse k-width of A5(k+q) is the size of the Cooper 
pair. This result is consistent with a general inequality 
derived in Ref. The authors of Ref. Q have fur- 
ther approximated the inequality to obtain A5(q) > 
(3/2)|Ap/2;, where A is the nearest-neighbor d-wave or- 
der parameter, and x is the concentration of holes. 

If Xm were of the magnitude predicted for a BCS su- 
perconductor, that is, iV(0)A, the intensity would be 
below the sensitivity limit of current neutron scattering 
measurements. That the intensity is readily detectable 
means that the condensate fraction in the high tempera- 
ture superconductors is more than an order of magnitude 
greater than that predicted in BCS theory. For tightly 
bound noninteracting Cooper pairs, acting as molecules, 
the condensate density would be the full density n. 

e. Underdoped superconductors: For underdoped su- 
perconductors Q, 5(k) is peaked above Tc at the 
wavevector q, but the growth of the intensity is gradual, 
increasing by about 20% from 300 K to Tc ^ 50 K, seen 
in both even and odd channels. In contrast, the intensity 
in the odd channel rises by almost a factor of 2 below 
Tc- It is very natural therefore to conclude that G(l) is 
of very different character than /(I). Thus, even in the 
underdoped case, it is possible to separate out the effect 
of G(l) and to focus on the contribution solely due to 
/(I). Before interpreting the data in terms of condensate 
fraction, we remark briefly on the peaked nature of iS(k) 
above Tc] there are a number of distinct possibilities. 

The first possibility is that the spin correlations are 
described by a quantum disordered state for which the 
structure factor is peaked at the commensurate wavevec- 
tor q. For such a quantum disordered state, it is difficult, 
however, to construct a d-wave modulation of the gap in 
the single particle excitation spectrum as seen in photoe- 
mission experiments |l0|] . 

The second possibility is superconducting fluctuations 
above T^. While such fluctuations must exist above T^,, 
especially in the underdoped regime, it is unlikely that 
they produce a peaked iS(k) at temperatures as high as 
300 K. 

The third possibility is a flux state for which the saddle 
point result [ p| for the insulator is remarkably isomor- 
phic to the BCS result for a superconductor, and iS(k) is 
peaked at the wavevector q, for essentially the same rea- 
son. Unfortunately, it is difficult to see why this should 
hold as the system is doped with holes, although it has 
been argued that the notions of the flux state should 
continue to hold for the doped case. 

The interpretation of A5(q), in terms of the conden- 
sate fraction has interesting consequences for the quan- 
tum phase transition as a function of the hole concentra- 
tion, X, hence for the global phase diagram. As x ^ x^ , 



typically ~5%, the superfluid density is known to vanish 
but the condensate fraction increases as shown in 
Fig. 1^. This is unusual and calls for an explanation. 

A possible resolution is that the experimental system 
cannot be considered pure, at least for x close to Xc, and 
that this transition is radically altered by disorder, allow- 
ing the condensate fraction to remain finite at the transi- 
tion. The discontinuous drop in the condensate fraction 
at Xc may suggest a first order transition, but it need 
not be. There is at least one example in statistical me- 
chanics where it is known that the order parameter van- 
ishes discontinuously, yet the correlation length diverges 
at the transition The role of disorder at this quan- 
tum phase transition may then be of crucial importance, 
and any discussion of a quantum critical point |^ must 
be reconsidered in this light. 
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FIG. 1. The imaginary part of the dynamical susceptibility, 
X+{q,i-^) of YBa2Cu306+i at 10 K [9]. This is proportional 
to AiS(q), because, at such low temperatures, the difference 
between the two is negligible. The subscript -I- refers to the 
enhancement over the normal state close to Tc; the data is 
for the odd channel, as defined in the text. The parameter 
S ~ 0.42 corresponds to the crtical hole concentration Xc at 
which the superconductivity disappears at T = 0. 

There may be another possible resolution The 
pseudogap phase may be a discrete broken symmetry 
phase in disguise, known variously as "fiux" or "dimer" 
ordered state, which has a peaked <S(k). Superconduc- 
tivity then develops on top of this. Such a broken sym- 
metry state is easily corrupted by disorder, notoriously 
plentiful in this class of materials, which is why it has 
been identified experimentally only as a "crossover" in 
some materials and not at all in others. This crossover 
will sharpen to an actual phase transition as the mate- 
rial qualities improve. The neutron intensity within the 
resonance peak is then a superposition of two effects, one 
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due to the condensate fraction and the other the due to 
the discrete broken symmetry. 

We cannot, of course, rule out the possibiUty that fu- 
ture measurements may reveal that the condensate frac- 
tion turns around and vanishes continuously at 
pecially because it is so difficult to obtain the absolute 
magnitude of the intensity in the neutron measurements. 
If this is the case, the interpretation of the quantum crit- 
ical point at Xc will be straightforward. There is some 
evidence that the intensity may vanish at Xc- The data 
shown in Fig. |l| was obtained by integrating over only the 
positive part of the imaginary part of the susceptibility 
enhanced over the normal state. Thus, a more complete 
integration over the frequency will be necessary to obtain 
the true theoretical value of the condensate fraction. 

/. Three-dimensional (3D) and bilayer couplings In 
bilayer superconductors, the adjacent layers are grouped 
into CuO-bilayers, which are then very weakly coupled. 
The small coupling between the bilayers should have a 
small effect on the 2D spin fluctuations seen above Tc, 
and indeed this is consistent with experiments. Likewise, 
such small couplings should have small effects also on the 
ordered state at T = 0. Close to Tc, small ?>D couplings 
should affect the critical behavior, however, and, below 
Tc, these couplings are necessary to establish an order 
parameter. As far as these weak ZD couplings of groups 
of bilayers are concerned, qualitatively, little needs to be 
added to the previous discussion. These comments would 
be hardly modified if we were considering superconduc- 
tors in which single layers (instead of groups of bilayers) 
were coupled by weak ZD coupling. 

In contrast, there is evidence that the individual layers 
within a bilayer are reasonably strongly coupled, thereby 
splitting the bilayer spin responses that are odd and 
even with respect to the individual layers. Although 
the dynamics is beyond the scope of the present pa- 
per, we can make a few brief qualitative remarks. For 
magnetic neutron measurements, the relevant coupling 
must be an antiferromagnetic superexchange, which must 
gap the even response function at low frequencies, as 
in experiments. This can be seen by noting that the 
spin correlations are approximately described by a bi- 
layer cr-model l]!^ and expanding the Euclidean La- 
grangian to quadratic order [p^ . The magnitude of 
the bilayer gap is a function of the cr-model param- 
eters and must depend on doping, as in experiments 
P|. As a result, from the ZD dynamic structure factor, 

the symmetric (even) part can be dropped at low frequen- 
cies. The integration across the resonance, which appears 
sharply in the odd (antisymmetric) channel, then gives 
the required structure factor. 

In conclusion, we have investigated the static struc- 
ture factor for which we could make precise statements, 
although the experiments on this are far from complete. 
We hope that future experiments will shed further light 



on the questions regarding the condensate fraction raised 
here. It remains to be seen if similar precise statements 
could be made for the dynamical structure factor as well. 
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